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> Abstract 

XT) 

ps) We perforin a phenomenological fit to all ATLAS, CMS, CDF and DO Higgs boson 

"^ data available after Moriond 2012. We allow all Higgs boson branching fractions, 

^ its couplings to standard model particles, as well as to an hypothetical invisible 

CN sector to vary freely, and determine their current favourite values. The standard 

model Higgs boson with a mass 125 GeV correctly predicts the average observed 

• wH rate and provides an acceptable global fit to data. However, better fits are obtained 

rN by non-standard scenarios that reproduce anomalies in the present data (more 77 

^ and less WW signals than expected) such as modified rates of loop processes or 

partial fermiophobia. We find that present data disfavours Higgs boson invisible 

decays. We consider implications for the standard model, for supersymmetric and 

fermiophobic Higgs bosons, for dark matter models, for warped extra-dimensions. 



1 Introduction 

Identifying the mechanism of electroweak symmetry breaking is the main goal of the Large 
Hadron Colhder (LHC). In the standard model (SM) the electroweak symmetry is broken due 
to the existence of an elementary scalar particle — the Higgs boson [1, 2, 3, 4, 5]. Based on 
data collected in 2011, both the ATLAS and CMS experiments at the LHC published results of 
their searches for the SM-like Higgs boson that, yet inconclusively, support its existence with 
a mass mh ^ 125 GeV [6, 7, 10]. Those results have been recently updated at the Moriond 
2012 conference, where all the Tevatron and LHC collaborations presented their updated Higgs 
boson searches as well as some new results. The combined Tevatron analysis of all collected 
data confirms the LHC excess around 125 GeV in the h —^ bb channel at 2.6cr level; CMS 
presented an improved 77 analysis; ATLAS presented new WW*, bb and rf searches with 
full 2011 luminosity. Furthermore, both the ATLAS and CMS experiments showed results of 
searches for a fermiophobic (FP) Higgs boson in the /i — )■ 77 channel that both show a positive 
hint around 125 GeV with local significances about 3a. This excess is consistent with the total 
inclusive 77 rate observed by the LHC [11]. 

Accidentally, nih ~ 125 GeV is a particularly fortunate value for the LHC, because, accord- 
ing to the SM predictions, various Higgs boson search channels are measurable. Those arise 
from a combination of SM Higgs boson branching fractions [8] 

BR{h -^ T+T-) = 6.4%, 
BR{h -^ 77) = 0.22%, (1) 

and production mechanisms with cross sections [9] 

a{pp -^h) = (15.3 ± 2.6) pb, a{jpp -^ jjh) = 1.2 pb, ,^. 

a{pp^Wh) = 0.57 ph, a{pp ^ Zh) = 0.32ph, ^' 

named gluon-gluon fusion {gg — )■ h), vector-boson fusion (VBF) and associated production 
with W and Z bosons (Vh). Because different search categories are sensitive to different Higgs 
boson couplings, the LHC can study the properties of a Higgs boson with nih ~ 125 GeV and 
test if it follows the SM predictions or is affected by new physics. 

With the presently collected statistics none of the search channels alone is sensitive to 
the SM Higgs boson nor are the combined results of Tevatron, ATLAS and CMS statistically 
conclusive. Therefore one expects large statistical fluctuations of the expected signal in all 
the search channels. Indeed, all measured LHC 77 rates, dominated by the new results in 
the VBF category, have central values above the SM prediction while all the WW* rates have 
central values consistently below the SM prediction. On the one hand, those anomalies may be 
statistical fluctuations. On the other hand, they may signal new physics beyond the SM. From 
a theoretical point of view, reconstructing the Higgs boson properties is an important way to 
address the main issue that LHC can clarify: is there a natural reason behind the the smallness 
of the weak scale, ruh <^ Mpi? Indeed, if the weak scale is naturally small, one expects that 
the new physics that cuts off the top loop contribution to ml (such as light stops at the weak 
scale in supersymmetric models) also affects the gg ^ h and h -^ 'j'j rates. Therefore a global 
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-» ZZ*) = 2.7%, 


BR(/i -^ gg) = 8.5%, 
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study of all the Higgs boson collider data obtained so far is necessary to test the SM and to 
discriminate between different new physics scenarios in the Higgs sector. 

In this work we study the collider data collected so far in Tevatron and the LHC in order 
to derive Higgs boson properties. The study of the Higgs boson fit was pioneered in Ref. [12], 
while recent fits of new LHC data were published in Refs. [13, 14, 15]. We improve on previous 
fits by including the new data presented in the Moriond 2012 conference, and by performing 
more general fits that cover a wider spectrum of new physics models. To achieve this goal 
we allow all the Higgs boson couplings to deviate independently from their SM values. We 
also allow for an additional Higgs boson invisible width, possibly due to decays into the dark 
matter. Anomalous features are dominated by the new results presented in Moriond 2012, 
disfavouring the SM compared to the previous fits and motivating new physics scenarios. We 
discuss implications of our results in the context of different models. More LHC data is needed 
to discriminate between those scenarios. 

The paper is organised as follows. In section 2 we describe the existing experimental results 
and the statistical procedure we adopt. In section 3 we perform the fits to data. In section 4 
we discuss implications of our results on different models. We conclude in section 5. 

2 Data and statistical analysis 

The experimental collaborations measure rates of Higgs boson signals R. Their results could 
be fully encoded in a likelihood C{R, rrih), but only a limited amount of information is reported 
by the experiments. All collaborations report the upper bounds on rates at 95% C.L., -Robserved, 
and the expected upper bound at 95% C.L. in absence of a Higgs boson signal, -Rexpected, as 
function of the Higgs boson mass rrih. Given that information, our statistical analyses follows 
the one outlined in Ref. [16]. Assuming that the x^ = — 21n£ has a Gaussian form in R, 

X^ = (i? - /i) Va^ (3) 

these two experimental informations allow one to extract the mean fi and the standard deviation 

fJ' -'''Observed -'•'expected) (^ Z 777^ ; \^) 

1.9b 
where 1.96 arises because 95% confidence level corresponds to about 2 standard deviations. 
The Gaussian approximation by construction agrees with the full result at this value of R, but 
away from it the approximation may be not accurate for channels that presently have a low 
number of events (such as /i — ?■ ZZ* — )■ Ai). We are aware of this fact, but at the moment it 
is difficult to do better using the available data. We verified that our procedure gives similar 
results as the refined procedure in [13] and that our procedure agrees better with present values 
oi fi± a, when reported by experiments at ruh = 125 GeV. 

We also neglect correlations of uncertainties among different measurements (e.g. uncertain- 
ties on luminosity and on the SM prediction, at the ±15% level and therefore subdominant 
with respect to present experimental uncertainties) and approximate the full x^ with 

x'=y: ^-^^^, (5) 



where the sum runs over all measured Higgs boson rates i. In the present stage of experimental 
accuracy such a simplified statistical framework captures the main features in data and allows 
us to study general properties of the data, that is the purpose of this work. 

We consider all available Higgs boson data reported at the Moriond 2012 conference and 
before: 

1. The pp — > Vh —7- Vhh rate measured by CDF, DO [17] and the related pp — )■ Vh — )■ Vhh 
rate measured by CMS and ATLAS [18]. 

2. The pp — )■ jjh — )■ jjWW rate measured by CMS [19]. 

3. The h -^ WW -^ 2£2z/ rates measured by CMS and ATLAS [20]. 

4. The h-> ZZ ^M rates measured by CMS and ATLAS [21]. 

5. The /i ^ 77 rates measured by ATLAS and CMS [22], and CDF, DO [23]. 

6. In the context of fermiophobic Higgs boson searches, CMS measured the pp — )■ jjh — )■ 
jJ77 rate [24] where the jj tagging is added to select Higgs boson produced via the 
VBF process. Indeed, the cuts performed by CMS {rrijj > 350 GeV, pxj^ > 30 GeV, 
Ptj2 > 20 GeV) significantly reduce the gg ^ h contribution, such that we estimate that 
reinterpreting this experimental result in a general context, it is roughly a measurement 
of 

[0.033a{pp ^h)+ a{pp -^ jjh)] x BR{h -^ 77) (6) 



with the result 



observed rate ^ ^ 

—^^-T- = 3.3 ±1.1 for m^ = 125 GeV. 7 

SM rate 



7. In the context of fermiophobic Higgs boson searches, ATLAS measured the pp — )■ hX — )■ 
77X rate with a high cut pTh > 40 GeV on the Higgs boson transverse momentum 
[25] (we are oversimplifying by omitting several secondary issues). This cut allows to 
suppress the gg ^ h production process, while keeping most of the signal in the VBF and 
associate production mechanisms. To see how much gg ^ h is suppressed we allowed for 
additional QCD jets performing simulations with the Pythia [26] and MadGraph [27] 
codes. We find that this experimental result can be re-interpreted in a general context as 
a measurement of 

[0.3a{pp ^h) + a{pp -^ Wh, Zhjjh)] x BR{h -^ 77), (8) 

with the result 

observed rate ^ ^ , ^ , , 

c,, , = 3.3 ±1.1 for m;, = 125 GeV. 9 

SM rate 

8. The h ^ TT rate as measured by CMS and ATLAS [28]. 

In the left panel of Fig. 1 we show our approximated combination of all Higgs boson data. 
Higgs boson masses around 125 GeV are favoured by the rate, and some ZZ and 77 events 
(which have little statistical power in fixing the rates but large resolution in mh) favor mh = 
125 GeV, the value that we will adopt in the rest of the paper. 
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Figure 1: Left; The Higgs boson rate favoured at la (dark blue) and 2a (light blue) in a global 
SM fit as function of the Higgs boson mass. Right; assuming m^ = 125 GeV, we show the 
measured Higgs boson rates at ATLAS, CMS, CDF, DO and their average (horizontal gray band 
at ±la). Here (red line) corresponds to no Higgs boson, 1 (green line) to the SM Higgs boson. 



Assuming rrih = 125 GeV, we summarise all data in the right panel of Fig. 1 together with 
their la error-bars, as derived by us following the above-described statistical procedure. The 
horizontal green line in the right panel of Fig. 1 is the SM prediction, and the horizontal red line 
is the background-only rate expected in the absence of a Higgs boson. The grey band shows the 
±1(T range for the weighted average of all data. It lies along the SM prediction. Furthermore, 
the global x^ of the SM fit is 17 for 15 dof. 

However, it is interesting to split data into three categories according to the final states and 
compute the average for each one of them: 



observed rate 
SM rate 



2.1 ± 0.5 photons 

0.5 ± 0.3 vectors: W and Z 

1.3 ± 0.5 fermions: h and r 



(10) 



This shows the main anomalous features in current measurements. First, the 77 channels 
exhibit some excess, mainly driven by the vector boson fusion data presented at the Moriond 
2012 conference. Second, there is a deficit in the vector channels. Finally, the average rate of 
fermionic channels lies along the SM prediction; here the new Tevatron combination for h ^ bb 
plays an important role. 

If the Higgs boson mass is different from 125 GeV, unlike what indicated by excess in 77 
distributions around this value of the invariant mass, then the 77 rate in eq. (10) would be 
reduced with respect to what we assume in our fit. 
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Figure 2: Predictions for the Higgs boson rates in different scenarios: SM, free branch- 
ing ratios of loop processes, free couplings, radion, top-phobic and fermiophobic, defined via 
eqs. (13,14,19). 

3 Reconstructing the Higgs boson properties 

3.1 Reconstructing the Higgs boson branching fractions 

The Higgs boson observables that can be most easily affected by new physics contributions are 
those that occur at loop level, the h — )■ 77, h ^ gg and gg -^ h rates. Because the latter two 
are related via CP, we use a common notation h ■^ gg to indicate both of them simultaneously. 
Those loop level processes are particularly relevant for the LHC Higgs boson searches because 
77 is the cleanest final state, and because gg ^ h is the dominant Higgs boson production 
mechanism. The left panel of Fig. 3 shows, as yellow contours with solid borders, the la and 
2(j ranges of a global fit to these two quantities in units of their SM predictions. The best fit 

corresponds to 

BR(h o gg) BR(h ^ 77) , 

!^ ^^^ ^03 ^^ -!^^ ^ 4 fill 

BR{h ^ gg)sM ' ' BR{h ^ 77)sm ' ^ ' 

that shows a significant deviation from the SM prediction — the first number allows to best 
fit the reduced WW* rates, and the second number allows to fit the enhanced 77 rates, in 
agreement with eq. (10). The x^ of the global fit is significantly lower with respect to SM, 
decreasing from 17 (for 15 dof within the SM) to 5.2 (for 13 dof in this more general fit). The 
black thick line in Fig. 2 shows the best-fit predictions for the various measured rates, allowing 
to see how the fit is improved. 

The gray region with dotted contours in Fig. 3 show the fit obtained omitting the 77 rates 
with cuts dedicated to vector-boson-fusion production (items 6 and 7 in the list above). In the 
latter case the agreement with the SM is improved showing that the such data category plays 
an important role in the fit. 
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Figure 3: Left; fit for the Higgs boson branching fraction to photons and gluons. In yellow with 
continuous contour-lines: global fit. In gray with dotted contour-lines: the fermiophobic Higgs 
boson searches are excluded from the data-set. Red dashed curve: the possible effect of extra top 
partners, such as the stops. Right; fits for the invisible Higgs boson branching fraction, under 
different model assumptions, as explained in section 3.2. 

3.2 Reconstructing the Higgs boson invisible width 

New physics can easily give a large effect providing an extra invisible [29] Higgs boson decay 
channel, for example into dark matter particles [30, 32]. Alternatively, the effective operator 
\dij_H'^H\^ similarly has the effect of rescaling all rates by a common factor [33]. 

In the SM total Higgs boson width is predicted to be r(/i)sM ~ 4.0 MeV at rrih = 125 GeV, 
too small to be measured directly. 

It is well known that measuring the Higgs boson total width at the LHC requires additional 
assumptions [12] because the gluon final states cannot be measured over huge QCD background. 
Let us explain how present data can probe the Higgs boson width, without directly measuring 
it. The gluon fusion production rates are proportional to T{gg -^ h). In view of approximate 
CP invariance we can assume that T{gg -^ h) = r{h — )■ gg) and we collectively denote them as 
h -H- gg. Then, one partial decay width can be reconstructed by data. By performing a global 
fit to the Higgs boson branching ratios in the context of theories where the decay widths are 
related we can reconstruct the total Higgs boson width. Of course this is based on theoretical 
assumptions, but the result gets significantly different only in highly deviant models, e.g. in 
models where the Higgs boson predominantly decays into light quarks (a decay mode not probed 
by present data). 

In order to emphasise the mild model-dependence of this fitting procedure we perform three 
fits under different assumptions. We show our results [x^ as function of the invisible branching 
ratio) in the right panel of Fig. 3. 



i) First, we perform a global fit of all data assuming the SM plus an additional invisible 
decay width, obtaining 

BR{h -^ invisible) = -0.1 ± 0.23; (12) 

ii) Next, we weaken the theoretical assumptions: we keep the h ^ gg and the /i — )■ 77 rates 
as free parameters, and marginalise with respect to them (red dashed curve). We see that, 
even without assuming the SM prediction for h o gg, a (weakened) bound on the Higgs 
boson invisible width can still be derived from present data; the best fit value becomes 
positive, but again the preference is not statistically significant. 

iii) Finally, we repeat the fit in i), but ignoring the data for 77 from the vector boson fusion 
channels, obtaining a weaker bound (dotted curve). 

Adding an invisible Higgs boson width has the effect of suppressing all observed rates and, 
according to our fit, this is not favoured by present data. 

3.3 Reconstructing the Higgs boson couplings 

In this subsection we extract from data the Higgs boson couplings to vectors and fermions, 
in order to see if they agree with the SM predictions. Trying to be as general as possible in 
describing the Higgs boson couplings, we proceed phenomenologically extracting from data the 
following parameters: 

^"^ „SM' -^^ „SM' -^* „,SM' -^^ „,SM' -^T „SM' ^ > 

i/W "V yt Ub iiT 

where gw is the WWh coupling; gz is the ZZh coupling, yt the top Yukawa coupling, yh the 
bottom Yukawa coupling and y-r the tau Yukawa coupling. All models considered in this work 
and presented in Fig. 2 are defined via eq. (13). The SM corresponds io Ri = 1 for all the 
couplings. These parameter Ri have the following effects: 

• the partonic cross sections for gg —> h and for gg — > tih get rescaled by R^; 

• the partonic cross sections for qq — )■ qqh and for qq — )• Vh get rescaled by Ry] 

• the decay widths h — )■ VV* get rescaled by Ry where V = {W, Z}; 

• the decay widhts h ^ ff get rescaled by Rj where / = {6, r, . . .}; 

• the decay width h — )■ 77, arising from the interference of one-loop diagrams mediated by 
the top and by the W, gets rescaled by {1.28Rw - 0.28RtY for rrih = 125 GeV; 

• similarly the decay width h — )■ Z7 (not yet measured) gets rescaled by (1.05i?^ — 0.05-Rj)^. 

A simplifying case considered in previous analyses [13, 14, 15] is a common rescaling factor 
a for Higgs boson coupling to vectors and a common rescaling factor c for Higgs boson coupling 
to fermions: 

a = Rv = Rw ^ Rzi c = Rt = Rb = Rt- (-'-4) 

8 
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Figure 4: Left; fit of the Higgs boson couplings assuming common reseating factors a and c 
with respect to the SM prediction for vector bosons and fermions, respectively. Right; fit to 
the t-quark and to b-quark and r-lepton Yukawa couplings assuming the SM couplings to gauge 
bosons. The best fit presently lies somehow away from the SM prediction, indicated in the figures 
as 'SM'. The point marked as 'FP' is the fermiophobic case, and 'Of denotes the top-phobic 
case. Negative values of the top Yukawa coupling are preferred because lead of an enhancement 
ofh-^-f-f. 



We show in the left panel of Fig. 4 the resulting fit (continuous yellow contours). For comparison 
the dashed contours show the result obtained ignoring the 77JJ data from CMS and ATLAS, as 
is also done in Fig. 3. This allows to compare our results with the ones of previous analyses [13, 
14, 15] (although some other data has also been modified and added by experiments). Our 
results essentially agree, up to the difference due to our use of more recent data. 

We see that a negative RtRw < is favoured because it implies a constructive interference 
between the top quark and W boson loops in the decays /i — )■ 77 increasing the corresponding 
rates. Notice that the new data prefers suppression of the WW* rates via suppression of the 
gg ^ h cross section, while the Higgs boson coupling to vectors can be somewhat larger than 
without VBF data. Notice also that the SM point (1,1) is disfavoured beyond 2a. 

In the right panel of Fig. 4 we assume the SM values for the Higgs boson gauge couplings 
[Rw = Rz = 1) and present a fit to the Yukawa couplings Rt and Rb = Rr- We, again, see 
that -Ri < if somehow favoured and the SM is disfavoured. The two branches approach the 
pure fermiophobic point (0,0), denoted by FP in Fig. 4, but pure fermiophobia is disfavoured 
by the fit. 

In Fig. 5 we consider the most general case where we allow all 4 parameters Rw = 
Rz, Rt, Rb, Rt to vary and show the favoured regions for the pairs Ry, Rt (left) and Rb, Rr 
(right) marginalised over the remaining two parameters. The main features of this global fit 
remain the same as in previous cases: RtRw < is favoured and Rw, Rb and R^ are constrained 
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Figure 5: Global fit for the Higgs boson couplings to vectors, to the t-quark, to the h-quark, to 
the T lepton. All these couplings are freely varied and in each panel we show the y^ as function 
of the parameters indicated on the axes, marginalised with respect to all other parameters. We 
again assume rrih = 125 GeV and find that the best fit presently lies somehow away from the 
SM prediction, indicated in the figures as 'SM'. 

to be around their SM values of 1. Fig. 2 shows the best fits (red dashed lines), both allowing 
for negative Yukawas (thick line) and restricting all Yukawas to be positive, as in the SM (thin 
line) . 

4 Implications for Higgs boson models 

In order to interpret our general results presented in Figs. 1-5 in the context of any particular 
model of Higgs boson, two logical possibilities arise. First, all the present anomalies in data, 
listed in eq. (10), could be just statistical fiuctuations. Second, the emerging pattern in eq. (10) 
could be real and signal new physics beyond the SM in the Higgs sector. Intermediate possi- 
bilities are of course possible. In order to discriminate between these possibilities, we present 
in Fig. 2 the predictions channel by channel of some particular scenarios that we studied in 
Figs. 3-5 for collider searches. The best fit x^ of those scenarios is also presented in the figure 
in order to compare different scenarios with each other. 

4.1 The Standard Model 

Naturally, the reference model for all comparisons in the previous sections is the SM. After 
fixing the Higgs boson mass to the best fit value rrih = 125 GeV, the SM does not have any 



10 



free parameter left to vary. Therefore all the anomalies in the present data must be statistical 
fluctuations and disappear with more statistics. This interpretation is supported by the fact 
that the average of all data agrees with the SM prediction, as seen in Fig. 2, and the global x^ 
is good: 17 for 15 dof (we recall that with n^ 1 degrees of freedom one expects x^ = ''^± ^A^)• 

On the other hand, our best fit (black curve in Fig. 2) has a significantly lower x^ = 5.5 
for 13 dof: a bigger reduction than what is typically obtained by adding two extra parameters 
(one expects Ax^ = —An ± vAn when adding An ^ 1 parameters). The SM is disfavoured at 
more than 95% CL in this particular context, but of course we added the two parameters that 
allow to fit the two most apparent anomalies in the data, the 77 excess and the WW* deficit. 

Only more data will tell if this is a trend, or if we are just fitting a statistical fluctuation. 

4.2 Fermiophobia and dysfermiophilia 

Fig. 2 shows predictions for different fermiophobic [34, 35, 11] scenarios. While bottomphobic 
Higgs boson is excluded by our fits, top-phobic or pure fermiophobic Higgs boson (with exactly 
vanishing Yukawa couplings) provide acceptable fits, of quality almost as good as the SM fit, 
despite that their predictions are significantly different. The pure fermiophobic model captures 
the features of data qualitatively correctly but predicts larger signal rates than is observed in 
the LHC, especially in the WW* + jj channel. In addition, the fermiophobic fit suffers from 
the h ^ bb signal claimed by Tevatron and CMS. 

The agreement of the fermiophobic Higgs boson with data can be improved by allowing a 
moderate small additional Higgs boson branching fraction, because this allows to reduce the too 
large prediction for the 77J j rate [1 1] , which is very sensitive to the precise value of the Higgs 
boson mass and width. In fermiophobic models such small Yukawa couplings can be generated 
via quantum effects [35] . We note that fermiophobia lowers the vacuum stability bound on the 
Higgs boson mass, allowing 125 GeV Higgs boson to be consistent with no new physics below 
Planck scale. 

Our fits in Figs. 4-5 show that reducing some or all of the SM Yukawa couplings allows 
to again significantly improve the global fit compared to the SM, down to x^ ~ 7. The main 
feature of the improved fit is yt ~ —0.7yf^, because this allows to enhance the /i — )■ 77 rate 
and reduce the gg -H- h rate. Admittedly, a 'wrong' Yukawa coupling to the top and to the 
other fermions (dysfermiophilia) is an even more serious pathology than fermiophobia. 

Various theoretical frameworks easily lead to modified Higgs boson couplings at moderate 
level. In models with more than one Higgs multiplet the Yukawa couplings of the light Higgs 
boson can be non-standard [36] (this is what can happen also in supersymmetric models). Al- 
ternatively, in models where SM fermions mix with extra fermions at the weak scale, integrating 
out the extra fermions, their effects get encoded in effective operators of the form ff'HH'^H, 
that lead to modified Higgs boson couphngs to the SM fermions / [37]. Such operators also 
arise in composite Higgs boson models. 

In models where the Higgs boson is a composite particle one generically expects that Higgs 
boson couplings get modified by form factors, approximated at low energy by effective opera- 
tors [33]. This is the framework considered in the fits of Refs. [13, 15]. In this kind of models, 
the rescaled SM expressions for these rates that we assumed remain valid even when new physics 
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is so large [41]. The Higgs boson coupling to W, Z vectors can be easily reduced by mixing the 
Higgs boson with other scalars; a good fit to electroweak precision data then demands that the 
extra scalars are not much heavier than the Higgs boson. Increasing the Higgs boson gauge 
couplings is theoretically more challenging [42]. 

4.3 Supersymmetry 

Supersymmetric theories that attempt to solve the naturalness problem of the electroweak scale 
have been stringently constrained by the LHC direct searches as well as by the Higgs boson 
results [44, 45]. Within the MSSM one needs light and strongly mixed stops, and there are two 
main modifications of Higgs physics. 

First, light stops modify the predictions and for the /i — )■ 77 and gg ^ h rates [38]. In 
practice their extra loop effect is described by a deviation of our parameter Rt from one: 



m^ 



R. = l+ ^ 



+ 



{At -/i/tan/3)^ 



2 ' 2 22 



(15) 



at leading order in the limit mi^ ^ '^ nit [14]. We see that the sign of the new effect is not fixed 
and can be negative in the presence of strong stop mixing. 

The red dashed curve in the left panel of Fig. 3 shows how these rates are affected by Rt 
(this applies not only to stops, but also to any extra particle with same gauge quantum numbers 
as the top, such as heavy top partners in little-Higgs models). Rt = 1 corresponds to the SM, 
and i?f = to the total suppression oi gg -k^ h. A Rt < 1 increases the /?,—)■ 77 rate, but 
only mildly because this rate is dominated by the W loop. The red dashed curve enters in the 
best-fit region when Rt ~ —1-7, a situation that cannot be achieved in view of bounds on the 
stop mass. We recall that such bounds are extremely model-dependent, because the signature 
depends on the unknown stop decay modes, and the production depends on the unknown gluino 
mass. For example, in gauge-mediated SUSY breaking the lower bound is 310 GeV [39]. If we 
assume m^^ ^ rrig, the bound on the stop mass is around 900 GeV [40], assuming that gluino 
decays always via a sbottom b into bNi with a neutralino mass ttitvi < 150 — 300 GeV. The 
bound on rrit gets about 200 GeV lower if the gluino decays fully via a stop. 

Second, due to the presence of two Higgs doublets Hi and H2, one has modification of the 
Higgs boson couplings at tree level. Our R parameters get modified as: 

sin ex cos cy 

Rw = Rz = sin(/3 - a), Rb = Rr = ^, Rt = ^—5-, (16) 

cos p sm p 

where tan/3 is the usual ratio between the two Higgs boson vev, and the a is the usual angle 
that diagonalises the mass matrix of Re {H^,!!^), with a — )■ /3 — 7r/2, in the decoupling limit. 
The angles a and P depend on the model and specific deviations arise depending on how 
ruh ~ 125 GeV is reached: extra D-terms imply an increase in /i — )■ 66 while extra F-terms lead 
to a decrease (unless extra singlets are light) [43]. The total Rt is the combination of the two 
effects discussed above. 

As previously discussed, both the WW* and the 77 rates are roughly proportional to R^; 
thereby this correlation prevents to go in the direction favoured by data (lower WW* and higher 
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77)5 cLS already observed in the context of numerical MSSM scans [44], and in extensions of the 
MSSM [45]. 

4.4 Dark matter models 

The main motivation for an invisible Higgs boson decay width comes from the existence of 
Dark Matter (DM) of the Universe. The Higgs portal [31] offers a natural possibility to couple 
the Higgs sector to the dark matter. If the dark matter particles are two times lighter than 
the Higgs boson, they can lead to invisible Higgs boson width. Because Higgs boson decays 
to fermion dark matter are essentially ruled out by direct detection constraints [32], in this 
scenario the dark matter is naturally scalar. 

Let us consider, for example, the simplest DM model obtained adding to the SM a real 
singlet scalar field 5* coupled to the Higgs doublet H by the —\S'^\H\'^ Lagrangian term [46]. 
Fixing the DM/Higgs boson coupling A assuming that the thermal relic DM abundance is equal 
to its cosmologically measured value I^dm = 0.112 ± 0.0056 [47] allows us to predict the Higgs 
boson invisible decay width and the direct DM detection cross section o"si 

V{h ^SS) = ^^ !l l^P'^, asi = ^'^'^^' . (17) 

8TTmh\j ml ' ^^ TrMl^ml ' 

The bound BRmv < 0.4 at 95% C.L. derived in section 3.2 then implies Mdm > 50 GeV and 
asi < 0.4 10~^"^ cm^, assuming the nucleon matrix element / = 0.3. While Higgs boson invisible 
decays to fermionic dark matter are already disfavoured, our work shows that also light scalar 
dark matter is not supported by data. 

4.5 Higgs boson or radion? 

The Higgs boson couples to the SM fermions with a strength proportional to fermion masses. 
Similar couplings can be obtained by considering an hypothetical particle (p, a radion, with a 
coupling to the trace of the SM energy-momentum tensor suppressed by some scale A: 

X^' = X (E ^/// - ^zZl - 2Ml,Wl + a\ . (18) 

In our language this is described by setting 

R = R^ =Rz = R^ = Rf^ = R^ = ^v/A, (19) 

where v = 174 GeV. One important difference arises at quantum level: scale invariance is 
anomalous, such that a new term appears. 






^ = '^^Gl^G^^ + -^F,.F,„ (20) 



where the numerical coefficients are the SM /3-function coefficients for the strong and electro- 
magnetic gauge couplings. As a result, the ip decay widths into gg and 77 differ from the 
corresponding Higgs boson decay widths [48]. 
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Such a particle is often called 'radion' because it arises in the context of models with one 
warped extra dimension as the mode that controls its size. However, in this kind of models the 
radion can appear with extra couplings and together with other unseen particles. 

We focus on the effective coupling in eq. (18), and find that the excess seen around 125 
GeV could be due to such a 'radion' rather than to the Higgs boson. The best fit is obtained 
at -R = 0.28 ± 0.03 (i.e. A ^ 870 GeV) and its quality is slightly worse that the best Higgs fit, 
as illustrated in Fig. 2. More data are needed to discriminate among the two possibilities. 

5 Conclusions 

We performed a global phenomenological analysis of all Higgs boson collider data available after 
Moriond 2012 (including those presented in the context of fermiophobic Higgs boson searches) 
assuming that the hints observed at rrih ~ 125 GeV arise from the Higgs boson. The SM 
provides an acceptable fit, however it is not favoured: present data with large uncertainties 
favor a /i — > 77 rate enhanced by a factor of ~ 4 and a gg ^ h rate reduced by a factor of 0.3. 
An invisible Higgs boson branching ratio larger than 0.4 is disfavoured, putting constraints on 
models where dark matter couples to the Higgs boson. Pure fermiophobic Higgs boson scenario 
gives almost as good fit as the SM but with significantly different predictions for the Higgs 
boson phenomenology. Partially fermiophobic scenarios are among those giving the best global 
fit. We find that the apparent excess can alternatively be interpreted as a 'radion' i.e. a particle 
similar to the Higgs boson, but coupled to the trace of the SM energy momentum tensor. 

More LHC data should clarify whether the present anomalies in data are statistical fluctu- 
ations or first evidence of physics beyond the SM. 
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